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We set up recursion relations for calculating all even moments of the end-to-end distance of a 
Porod-Kratky wormlike chains in D dimensions. From these moments we derive a simple analytic 
expression for the end-to-end distribution in three dimensions valid for all peristence lengths. It is in 
excellent agreement with Monte Carlo data for stiff chains and goes properly over into the Gaussian 
random-walk distributions for low stiffness. 



I. INTRODUCTION 

The main quantity of interest for a stiff polymer in D 
dimensions is the end-to-end distribution defined by the 
path integral of a stiff polymer |]| 

P L (R)cx J d D u b J d D u a J V D u 

X8W (li- £ dsn(s)j e -t*/*)J L *>[*'(>)]'. (1.1) 

where R is the reduced stiffness related to the persistence 
length £ by 



k R T 



(1.2) 



The unit vectors u(s) are the tangent vectors of the space 
curve of the polymer. A Fourier representation of the 5- 
function brings this to the form 

pica jD\ _ f f 

P L (R)oc/ ^e hXR / 2 J d D u b J d D u a (u b L\u a 0) x ,(l.3) 



where 
(u & L|u a 0) 



u(L)— Uf, 

A / D D ue -(*/ 2 )/» t ' is {i u '( s )] ! + A - u ( s )} 



(0)=u Q 



(1.4) 



coincides with the euclidean path integral of a point par- 
ticle of mass M = R moving on a unit sphere in an ex- 
ternal electric field A. 

Since u(s) are unit vectors, the path integral is not 
solvable exactly, except for zero A. It is however, easy 
to find arbitrarily high even moments of the end-to-end 
distance of the distribution Pj y (R): 



(R 2n ) = J d D RR 2n P L (K). 



(1.5) 



The 2nth moment of the chain can be obtained directly 
from the expansion coefficient in powers of A of the inte- 
gral over Q1.4|) . 



In naturla units with «)1, the path integral ljl.4|l solves 
the Schrodinger equation 



-iA u + iA.u + £)mru.,Or A = 0. il.(is 



where A is the Laplacian on a unit sphere. The exter- 
nal electric field A may be assumed to point in the z- 
direction, or the Dth direction in D-dimensions. In the 
distribution (|1.3|l . only the integrated expression 



i>{z,T)\) = / d D u a (ur|u a 0) ; 



(1.7) 



appears, which is a function of z = cos 9 only, where 6 
is the angle between u and the electric field A. Setting 
R = 1, for a moment, we obtain the simple Schrodinger 
equation in euclidean time 



Hip{z, t; A) = -—Tp{z, r; A), 



(1.8) 



H = H + A Hj 



1 A 1, 



^-z 2 )^-{D~l)z- 



dz 



dz 



+ -\z (1.9) 



Now the desired moments 111 ,')H can be obtained from 
the coefficient of X 2n /2 2n (2n)\ in the expansion of the 
integral over i|1.7|) . 



/(t; A) = / dziJj(z,T;X), 



(1.10) 



in powers of A, evaluated at the euclidean time t = L. 



II. RECURSIVE SOLUTION OF THE 
SCHROEDINGER EQUATION. 

The function /(£<; A) has a spectral representation 

~ i^VOt( )ex P (-EWL) g(0(0) 
J{L; A) = y j — ; 77777^ 7777— A ZA ) 
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where ip^(z) are the solutions of the time-independent 
Schrodinger equation H<p^'(z) = E^ tp® (z) . Apply- 
ing perturbation theory to this problem, we start from 
the eigenstates of the unperturbed Hamiltonian H = 
—A/2, which are given by the Gegenbauer polynomials 

Cf /2_1 (z) with the eigenvalues = l(l + D- 2)/2 0. 
Next we set up a recursion scheme for the perturba- 
tion expansion of the eigenvalues and eigenfunctions as 
described in We begin with a brief review of the 
method. Starting point is the usual expansion of energy 
eigenvalues and states in powers of the coupling constant 
A: 



3=0 

oo 

\r- ■)= E -rgA^'IO ■ 

V, i=0 



,(0\ 



(2.2) 
(2.3) 



The wave functions ip"> (z) are the scalar products 
(z\<p^). The index i counts the order of the interac- 
tion strength A. The lowest expansion coefficients of the 

energy are of course = Eq . In the second line, we 
have introduced auxiliary normalization constants ay for 
convenience to be fixed later. The state vectors |Z) of 
the unperturbed system are normalized to unity, but the 
state vectors \<*p^) of the interacting system will be nor- 
malized in such a way, that (<p {l) \l) = ai holds to all 
orders, implying that 



(0 
111 



= 5 



i',0 



(0 X 

7 fe ,o = °i,k 



(2.4) 



Inserting the above expansions into the Schrodinger 
equation, projecting the result onto the base vector (k\a.k, 
and extracting the coefficient of X 1 , we obtain the rela- 
tion: 



3=0 



Oik 



3=0 



where Vkj = X(k\z\j) are the matrix elemts of the inter- 
action between unperturbed states. For i = 0, Eq. (|2.5(l 
is satisfied identically. For i > 0, it leads to the following 
two recursion relations, one for k = I: 



(2.6) 



n=±l 



the other for k ^ I: 



'k,i (k) (I) ' 



where only n = — 1 and n = 1 contribute to the sums 
over n since 

(/) ^ai+n ^ + ^ = Q ^ for n ^ ±L (2 8) 



The vanishing of Wn^ for n 7^ ±1 is due to the band- 
diagonal form of the matrix of the interaction z in the 
unperturbed basis It is this property which makes 
the sums in 1)2. 6|) and (|2.7|) finite and leads to recursion 
relations with a finite number of terms for all e„- and ™^ 



To calculate Wn , it is convenient to express (l\ z \l + n) 
as matrix elements between unnormalized noninteracting 
states \l'} as 



(l\z\l + n) 



{l\z\l + n} 



^{l\l}{l + n\l + n} 



(2.9) 



where 



{k\F(z)\l}^ j dzC^-\z)F{z)C^-\z){l-z^^ 

(2.10) 

yielding [j| 

2 4 - D T(1 + D-2)tt 



l\ (21 + D -2) r(D/2-l) 2 



(2.11) 



Expanding the numerator of 1)2. 9|) with the help of the 
recursion relation for the Gegenbauer polynomials [|| 

(I + 1)|Z + 1} = {21 + D - 2) z \l} - (I + D - 3)|i - 1} 

(2.12) 

we find the only non-vanishing matrix elements to be 
l + l 



{l + l\z\l} 



21 + D -2 



{Z + l|/ + l}, 



(2.13) 



^~ 1 W- 2i + £? zl ^- 1 ^- 1 )' ( 2 - 14 ) 
Inserting these together with (|2.11|) into (|2.9|l gives 



(l\z\l-l) 



1(1 + D -3) 



(2Z + L>-2)(2Z + L>-4) 



(2.15) 



and a corresponding result for (Z|z|Z + 1). We now fix the 
normalization constants aj< by setting 



W (0 = 2!±l < ^|i + i ) = i 
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for all I, which determines the ratios 



07 



= (l\z\l + l) = 



1 (l + l) (l + D-2) 
ai+i • ■ ■ y (21 + D) (21 + D-2)' 

Setting further a,\ = 1, we obtain 

1/2 



(2.16) 



(2.17) 



a/ = 



n 

3=1 



(21 + D -2)(2l + D -A) 
1(1 + D-Z) 



(2.18) 
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Using this we find from (|2.8() the remaining nonzero Wn* 
for n = — 1: 



W 



1(1 + D -3) 



(2Z + D-2)(2Z + D-4) 



(2.19) 



We are now ready to solve the recursion relations (|2.f)() 
and (|2.7|l for jrl and order by order in i. For 

the initial order i = 0, the values of the 7^ are given 
by Eq. (|2.4|l . The coefficients ef^ are equal to the 



III. FROM MOMENTS TO END-TO-END 
DISTRIBUTION IN D = 3 DIMENSIONS. 



The moments can now be used to recover the exper- 
imentally accessible end-to-end distribution of the poly- 
mers for various degrees of stiffness. We parameterize 
the distribution with an analytic form 



Pl(R) ocr fc+2 (l-r p ) 



unperturbed energies = Etf' = 1(1 + D - 2)/2. For 



each 



1,2,3,..., there is only a finite number of whose moments are 



non-vanishing 7^ and e^- with j < i on the right-hand 

sides of (|2.6|l and l|2.7() which allows us to calculate yjp. 

and on the left-hand sides. In this way it is easy to 
find the perturbation expansions for the energy and the 
wave functions to high orders. 

Inserting the resulting expansions (|2.2|) and (|2.3|l into 
Eq. (|2.1|) . only the totally symmetric parts in ip^(z) will 
survive the integration in the numerators, i.e., we may 
insert only 



cz> W (z) 

rsymm J 



(^K ) mm )=E^SA i ^io>. 



(2.20) 



The denominators 

„« .12 \2 



of (|2.1() become explicitly 
i \lv'i a i'\ 2 ^ 2 *' where the summation over i is 
limited by power of A 2 up to which we want to carry 
the perturbation series; also I' is restricted to a finite 
number of terms only, because of the band-diagonal 

structure of the 7}/^. 

Extracting the coefficients of the power expansion in A 
from 1)2. If) we obtain all desired moments of the end-to- 
end distribution, the lowest two being, after reinserting 
R from (|1.2|) . 



{Ri) = 2{t;L-e 



1 



-L/i 



(2.21) 



(R 2n ) = 



= R/L. 



1) 



(3.1) 



r(-±*)ro 



i) 



(3.2) 



We now adjust the three parameters k, f3, and m to fit 
the three most important moments of the distribution 
exactly, ignoring all others. If the distances were dis- 
tributed uniformly over the interval r G [0, 1] , the mo- 
ments would be (i? 2 ™) = l/(2n + 2). Comparing 
our exact moments (i? 2 ™) (£) with the flat ones, we find 

that (i? 2n ) (£)/ (i? 2ll ) 9at has a maximum for n close to 
^max(C) = 4£- The most important moments turn out to 
be the ones with n = ?"i max (£) and n = n. max (£) ± 1. If 
(0 < 1 we use the lowest even moments (R 2 ), (R 4 ), 
and (i? 6 )- With these adjustments the resulting distribu- 
tions are shown in Fig. ^for various persistence lengths 
£. They are in excellent agreement with the Monte 
Carlo data (symbols) obtained by Wilhelm and Frey 0, 
and better than their one-loop perturbative results (thin 
curves) which are good only for very stiff polymers. For 
the small persistence lengths £ = 1/400, 1/100, 1/30, 
the curves are well approximated by Gaussian random 
chain distributions on a lattice with lattice constant 
a c ff = 2£ which ensures that a e ff = 2£ the lowest mo- 
ments (R 2 ) = a c f[L are properly fitted. 



<i? 4 ) 



4(£>+2) 



D 



D 2 + 6D-1 D-7 



D 2 



-L/i 



D+l 



4£ 4 



23D 2 -7D + 1 ^(D 



D 3 



(D + l) 2 



-L/i 



(D-5) 5 2DL/(D-l)i 

D*(D + iy 



The calculation of higher moments is straightforward 
with a Mathematica program, which we have made 
available on the internet in notebook form Ml . 
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FIG. 1: Distribution of the end-to-end distances of poly- 
mer for different stiffnesses, parametrized by the persistence 
lengths £ = 1/400, 1/100, 1/30, 1/10, 1/5, 1/2, 1, 2. They 
are compared with the Monte Carlo calculations of Wilhelm 
and Frey |J (symbols) and with his large-stiffness one-loop 
perturbative results (thin curves). For the small stiffnesses 
£ = 1/400, 1/100, 1/30, the curves are well approximated by 
Gaussian random chain distributions on a lattice with lattice 
constant a g = 2£ which ensures that a c g = 2£ the lowest 
moments (R 2 ) = a e ftL are properly fitted (dashed curves). 
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